/**
 * @file   initial_state.hpp
 * @author Kimmo Luoma <kimmo@PHOTON3>
 * @date   Fri Sep 14 13:33:59 2012
 * 
 * @brief  Class for the physical system definitions
 * 
 * 
 */
#ifndef PHYSICAL_SYSTEM_HPP
#define PHYSICAL_SYSTEM_HPP

#include<boost/numeric/ublas/vector.hpp>
#include<complex>
using namespace boost::numeric::ublas;

class initial_state{

private:
  vector<std::complex<double> > Psi0_;
public:
  initial_state();
  initial_state(vector<std::complex<double> > psi);
  initial_state(const initial_state& is);
  initial_state& operator=(const initial_state& is);
  ~initial_state();
  
  vector<std::complex<double> > operator()();
  
};

#endif
